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SUPERCRITICAL MINIMUM MEAN-WEIGHT CYCLES 


JIAN DING*, NIKE SUN^ AND DAVID B. WILSON 


Abstract. We study the weight and length of the minimum mean-weight cycle in the 
stochastic mean-field distance model, i.e., in the complete graph on n vertices with edges 
weighted by independent exponential random variables. Mathieu and Wilson showed that 
the minimum mean-weight cycle exhibits one of two distinct behaviors, according to whether 
its mean weight is smaller or larger than l/(ne); and that both scenarios occur with positive 
probability in the limit n —>■ co. If the mean weight is < l/(ne), the length is of constant 
order. If the mean weight is > l/{ne), it is concentrated just above l/(ne), and the length 
diverges with n. The analysis of Mathieu-Wilson gives a detailed characterization of the 
subcritical regime, including the (non-degenerate) limiting distributions of the weight and 
length, but leaves open the supercritical behavior. We determine the asymptotics for the 
supercritical regime, showing that with high probability, the minimum mean weight is 
(ne)“^[l + 7r^/(21og^n) -I- 0((logn)“^)], and the cycle achieving this minimum has length 
on the order of (logn)^. 


1. Introduction 

Given a directed or undirected graph with edge weights, a minimum mean-weight eycle 
(mmwc) is any cycle that minimizes the mean weight (ratio of total weight to cycle length) 
over all cycles in the graph. Finding an MMWC is a fundamental subproblem to a wide variety 
of algorithms: for example, they have been used ([GT89, RG94], building on [Kle67]) to give 
a strongly polynomial-time algorithm for the minimum-cost circulation problem, of which the 
maximum flow problem is a special case. Other applications of MMWC’s include algorithms for 
multicommodity flow [OMOO] and asymmetric travelling salesman tours [KW98]; for further 
applications see the extensive discussion in [Das04]. 

Several mmwc algorithms are available (see [Das04, GGTW09] and references therein); 
and it is of practical interest to understand their runtime in “average-case” settings, that is to 
say, on random inputs. Experimental studies [DG98, DIG99, Das04, GGTW09] suggest that 
an algorithm due to Young-Tarjan-Orlin [YT091] (based on an improvement of a parametric 
shortest-path algorithm [K081]) has the best runtime in standard random graph ensembles, 
where it substantially outperforms its worst-case theoretical guarantees. 

Motivated by the empirical studies, a natural direction is to understand the typical behavior 
of the MMWC in random graphs. With this in mind, Mathieu-Wilson [MW13] study the 
mmwc in the stochastic mean-held distance model: the complete graph, or complete digraph, 
with i.i.d. random edge weights. This is a canonical setting for the study of combinatorial 
optimization problems, both in the mathematics and physics literature: other examples include 
minimum spanning tree [Fri85, FM89], shortest path [FG85, Jan99, HHVM07, HVM08], 
traveling salesman [MP86, Fri04, WaslO], assignment [MP85, AldOl, LW04, NPS05], spanners 
[GFMS09], and Steiner tree [BGRS04, AFW12]. The analysis of [MW13] uncovers an unusual 
dichotomy for the MMWC in the stochastic mean-held distance model, and characterizes 
the subcritical regime. In this work we complement their analysis by characterizing the 
supercritical regime. 

Research supported in part by *NSF grant DMS-I313596 and ^NSF MSPRF grant DMS-I40II23. 
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1.1. Main result. Consider the stochastic mean-held distance model where each edge is 
independently weighted by an exponential random variable of unit rate. In this setting we 
study the minimum mean-weight cycle, which is unique with probability one. Denote its mean 
weight by and its length by so the cycle has total weight This model exhibits 

an unusual phase transition [MW13] at nWn = 1/e: In the subcritical regime nWn < 1/e, 
nWn is non-concentrated, and stays bounded; moreover the limiting distributions of 
nWn and are precisely characterized. In the supercritical regime nWn > 1/e, nWn is 
very concentrated with nWn — 1/e ^ 0 in probability, while diverges with n, at least 
on the order of (log log log n. The conclusions of [MW13] are much less precise in the 
supercritical case, leaving open the asymptotic order of nWn — 1/e and 

In this paper we characterize the asymptotics of and in the supercritical regime. 
In particular, our results conhrm simulations done by Uri Zwick and the last author which 
suggested that — (logn)^. Our main theorem is as follows: 


Theorem 1.1. In the complete graph or complete digraph with i.i.d. unit-rate exponential 
edge weights, let Wn and be the mean weight and length of the minimum mean-weight 
cycle, and write 


c„{n) = e ^[1-t-7r7[2(logn)^]]. 
For all e > 0 there exists a constant C = C{e) > 0 such that 

\n'Wn - c*(n)|(logn')^ < C 


lim inf P 

n^oo 


and 1/C < .if„/(logn)^ ^ C 


nWr, > e ^1^1 


e. 


( 1 ) 

( 2 ) 


These bounds are optimal in the sense that for any interval I with length |J| ^ I/O, 


lim sup P 

n^oo 


[nWn - c*(n)](logn)^ e I 
or Jd’n/i^ogn)^ e I 



( 3 ) 


Theorem 1.1 implies that in the supercritical regime nWn > 1/e, the random variables 
log[=Zfn/(logn)^] and (logn)^[n#k — c*] are tight, with non-degenerate distributions. 


1.2. Proof ideas. In the remainder of this introductory section we highlight some of the main 
new ideas in our proof, which allow us to overcome obstacles in the analysis of [MW13]. We 
remark that very similar obstacles arose in an analysis [Dinl3] for a related model, percolation 
of averages, which asks for the longest path with mean weight below a parameter A. While 
our work is carried out in the context of the MMWC, we expect our methods may be applied 
to improve results in [Dinl3]. 

In order to show that the supercritical MMWC has mean weight Wn in some interval [a, b] 
and satishes some property P, one must show (i) there are no cycles with mean weight in 
[l/(ne),a], (ii) there are no cycles violating P in [a, 6], and (iii) there is at least one cycle 
(satisfying p) in \a,h\. In particular, (ii) is usually done by hrst moment arguments, and is 
clearly easier for less restrictive properties P. On the other hand, a natural approach for (iii) 
is the second moment method, for which it is often advantageous to make P more restrictive. 

Indeed, it is demonstrated in [Dinl3, MW13] that a straightforward second moment method 
on the number of cycles fails, due to an excessive contribution from atypically light cycles (or 
paths) — conditioned on finding one atypically light cycle it is very likely to have a large 
number of light cycles overlapping with ^. To address this issue, these works developed a 
notion of uniformity: for a k-cycle with edge weights wi,... ,Wk summing to kc/n, we say the 
cycle is y4-uniform if the process X with increments W+i — Xi = {nwi/c— 1) has range at 
most A — meaning the cycle has no excessively light subpath. The typical MMWC is uniform 
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with high probability; but on the other hand the count of uniform cycles has small enough 
variance for the second moment method to go through. 

However, this uniformity property is not sufficiently restrictive to yield accurate implications 
on the asymptotics of Wni ^n- Indeed, it is clear that the precision achievable for the window 
[a, h\ is dictated by the accuracy with which P captures the typical properties of the MMWC. 
In this work, we restrict further to uniform cycles with a typical profile. We defer the formal 
dehnitions to Section 5; roughly speaking, we restrict to cycles for which the associated 
process X not only has range < A, but furthermore has typical local times within its range. 
The analysis of Sections 3 and 4 will show that this restricted property is satished with 
high probability. On the other hand, in Section 5 we show that the restriction captures the 
remaining variance, so that the number of cycles restricted in this manner is well-concentrated 
about its mean. Theorem 1.1 follows as a consequence. 

A key technical ingredient in our proof is a collection of precise estimates for exp-minus- 
one random walks (that is, a random walk with increments distributed as — 1, with E 
a unit-rate exponential variable) conditioned to have restricted range. While our results 
are reminiscent of analogous estimates for simple random walk or Brownian motion, for 
general random walks under suitable moment assumptions there is no general theory yielding 
estimates to the level of accuracy needed for our MMWC analysis. Our estimates for the 
exp-minus-one walk are derived in Section 3. A crucial input to these estimates is a precise 
characterization of the principal eigenvalue and eigenfunction for the exp-minus-one walk in 
an interval with absorbing boundaries. This analysis may be of independent interest, and is 
presented in Section 6. 

Acknowledgements. Computer experiments on the mean-held stochastic MMWC performed 
by Uri Zwick with D.B.W. provided valuable intuition at an early stage of this project. J.D. 
and N.S. thank the MSR Redmond Theory Group for its hospitality. 

2. Preliminaries 

Notation. We write fn ^ Qn (or pn ^ /n) if there exists an absolute constant G > 0 such 
that fn < Cpn for all n e N. We write ~ to indicate that < and > both hold. For numerous 
parameters in the paper, we write for example C = G(A, e) to indicate that G > 0 is a number 
depending only on A and e. 

2.f. Cycles and paths. For a cycle ^ we write len(‘^) for the length (number of edges) 
of when len('^) = k we refer to ^ as a fc-cycle. We write wgt('^) for its total weight, and 

c(^) = nwgt(^)/len(^) (4) 

for its mean weight scaled by n; we say ^ is c-light if c(^) ^ c. We apply these terms to 
paths as well as cycles; note that a fc-cycle involves k vertices while a fc-path involves k + 1 
vertices. To treat both undirected and directed random networks in a fairly unihed manner, 
we will always take cycles and paths to be directed. In our random network, the weight of 
any given fc-cycle or fc-path is distributed as the sum of k independent unit-rate exponential 
random variables: that is to say, a gamma random variable with shape parameter k, with 
probability density 

p—x^k—l 

= (fcTri)! ’ ^ ^ 
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We abbreviate this distribution as Gam(A;). For the sake of review, we repeat the following 
calculation from [MW 13]: 


Lemma 2.1. Let count all c-light k-cycles (k ^ 2), and let count all c-light k-paths 
(k ^ 1). For all c < 1 and 0 < 6 < 1 we have 





(n)t (ce)*[l - (1 - f)‘’] 

nk /^ 3/2 




{n)k+i {cef {cef 

-^ 77 - 

n^ /ci/2 /jl/2 





Proof. In the complete graph on n vertices, the number of (directed) /c-cycles (k ^ 2) is 
{n)k/k. The weight of any given fc-cycle is distributed as a Gam(fc) random variable, therefore 

= (!^p(Gam(fc) ^ ck/n) = f ^ 

k ^ fc Jo (k-iy. 

Since c < 1 we have — 1 uniformly over the range of integration, therefore 






rckjn 


X 


k-1 


dx 


{n)k {ck/n)^[l - (1 - 5)^] 


k 


k\ 


^ Jc{l—5)k/n !)• 

{n)k {ce)X (ce)^[l - (1 - 6)^] 


^3/2 ^3/2 ’ 

proving the estimate for c-light /c-cycles. The estimate for c-light fc-paths follows by noting 
that the number of fc-paths {k ^ 1) is (n)fc+i. □ 


An easy variation of the preceding calculation shows that in the targeted regime for 
there can be no cycles of length less than (logn)^: 


Lemma 2.2. Given any c = l/e[l -I- 0(l/(logn)^)] and G > 0, the probability that there is 
a cycle ^ with 1/e < c(^) ^ c and len(^) < G(logn)^ tends to zero in the limit n —» oo. 


Proof. The expected number of directed /c-cycles that are c-light but not (l/e)-light is 





{n)k e ^ ^ (ce)*^ — 1 ^ fc/(logn)^ 

^ k/ien) W^- ^ ^ ’ 


where the last bound used that k < (logn)^. Summing over 2 < fc < (logn)^ and applying 
Markov’s inequality proves the claim. □ 


2.2. Uniformity. It was previously demonstrated in [Dinl3, MW13] that cycles which are 
uniform (in a sense dehned formally below) have low variance. For a sequence of weights 
tci ,... ,Wk summing to kc/n, let us dehne its excedance relative to c (for short, c-excedance) 
to be the quantity 

^{nwi/c-l) = k{c/c-l). (6) 

i 

Definition 2.3. Say that a cycle ^ is (c, A)-uniform if it has no subpaths with c-excedance 
outside [—A, A]. We say simply that ^ is A-uniform if it is (c(^), A)-uniform. We apply 
the same terminology to paths as well as cycles. 
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Definition 2.4. Given weights wi, ... ,Wk with mean c/n, the D-tilted bridge is the process 
j 

Wj = — c — D/k), started from Wq = 0 and ending at W4 = —D. 

We refer to the D = 0 case as the untilted bridge: in particular, a cycle with mean weight c/n 
is A-uniform if and only if its untilted bridge has range < cA. Sometimes we say fc-bridge to 
emphasize that the bridge is defined on the time interval [0, fc]. 

2.3. Exponential random walks. The exp-minus-one random walk is the random walk on 
the real line with step distribution u — 1, where u is a unit-rate exponential random variable. 
Observe that the cycle bridges (Definition 2.4) can be rescaled to an exp-minus-one version: 
to see this, take any positive /i and consider the random vector w_ = (tci,... ,Wk) where the 
Wi are independent exponential random variables with mean /i. Following the notation (4), 
suppose they sum to kc/n. Let s be any positive c-measurable random variable, and consider 
the process 

i j 

Wj = — s) = s'^{nwi/s — 1), 0 < j < k. 

Conditioned on c, the vector w is distributed as a uniform sample from the space of all 
non-negative vectors in with sum kc/n, regardless of the value of fi. In particular, taking 
H = s/n shows that 

the process W/s is distributed as an exp-minus-one random walk . . 

started from the origin and conditioned to be at k{c/s — 1) at time k. ^ ^ 

For example, by taking s = c, we see that the process with increments {nwi/c — 1) is dis¬ 
tributed simply as an exp-minus-one walk conditioned to return to the origin at time k. We 
will apply (7) with s ^ c in the proof of Lemma 4.4. 

In view of Definition 2.4, we study the exp-minus-walk with restricted range. In Section 6 
we give a precise computation of the principal eigenvalue of the exp-minus-walk with 
killing outside [0, A]. In the limit of large A it behaves as 

XA = exp{-nyi2A^) + 0{l/A^)}. (8) 

By comparison, for simple symmetric random walk on the integers killed outside {1,..., A — 1} 
(with A integral), the principal eigenvalue is = cos(7r/74), which behaves in the limit of 
large A as exp{—7r^/(2y4^) -I- 0(1/^^)} (see e.g. [Kac45]). 

2.4. Proof overview. Having finished our preliminary calculations, we conclude this section 
by outlining the proof of our main result. 

In Section 3 we prove the necessary estimates for (range-restricted) exp-minus-one walks — 
taking as input the principal eigenvalue and eigenfunction for the walk with killing outside 
the interval [0,^4], which will be computed in Section 6. The most important consequences 
of this section concern the exp-minus-one random walk k-bridge, by which we mean an 
exp-minus-one random walk conditioned to return to the origin in k steps. Lemma 3.9 
computes the probability R\ for this process to have range < A: 

k3/2 

Ra ^ (Aa)"^ for all k > Hi 

Lemma 3.10 shows that if we condition this process to have range < A, then its local times 
are comparable with those of the analogously range-restricted Brownian motion. 
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In Section 4 we apply the random walk estimates to rule out supercritical cycles which are 
atypically light or long. Recall that in Lemma 2.1 we computed the expectation of the number 
of c-light fc-cycles. Let Z^{A) count c-light fc-cycles that are A-uniform: in Lemma 4.3 
we apply the estimate on R\ to prove 

EZ^{A) = iEZ^)R\ ^ for all k > A^. 

If a c-light cycle fails to be A-uniform, then we can extract a subpath whose bridge decreases 
by —A + 0(1) and is {A + 0(l))-uniform. We let Z^[A) count Opaths of this type; in 
Lemma 4.4 we apply random walk estimates from Section 3 to show that 


EZf(kl) < 


(ccAa)^ n 
A3 ^ 


for all i ^ 1. 


Consequently, if Zf^^ counts all c-light cycles of length ^ we have 

(ccAa)^ " 

k^l 


> 0) < 2 EZ^A) + Y,e.Z'M) S L 




n 

1 + ^ 
e^J 


( 9 ) 


This calculation suggests that we take A = logn -I- 0(1) and rule out values of c that make 
ccXa too small. This leads to the definitions 


Ao = log n, 



1 

e 


1 + 


2 (log n) 2 


+ 0(l/(logn)’) , 


( 10 ) 


applying (8). It suffices to prove Theorem 1.1 with Co in place of c*. The lower bound on 
nWn stated in the theorem is an easy consequence of Lemma 2.2 and (9), and this is the hrst 
main consequence of Section 4. 

The second main part of Section 4 is to rule out cycles that are much longer than (logn)^ 
in the regime |c — Co| l/(logn)3. This argument is rather more involved, but it is guided by 
the same basic principle that the untilted bridge of a cycle either stays in a restricted range, 
or has sharp decreases over restricted ranges. If a cycle is very long, its bridge must either 
have many sharp decreases, or else stay in restricted ranges over long intervals. Both events 
impose a severe probability cost which can be used to rule out the presence of long cycles. 

In Section 5 we identify a subcollection of uniform cycles with a “typical local time prohle,” 
and show that the number of such cycles is well concentrated about its mean. It follows 
that these cycles exist (with high probability) whenever their expected number is large. It 
follows from Lemma 3.10 that most of the contribution to EZ^(A) comes from cycles having 
a typical local time prohle — meaning that EZ^(A) locates the sharp transition. Theorem 1.1 
then follows in a straightforward manner. 


3. Random walk estimates 

In this section we derive the necessary estimates for exp-minus-one walks subject to 
restrictions on the range of the walk. A few of the estimates require some understanding of 
the principal eigenvalue Aa, and associated left eigenfunction (pA, for the walk with killing 
outside the interval [0, A]. These will be computed in Section 6; in the present section we 
shall require only the facts that 

Aa = exp{-[l -L 0(l/A)]7r7(2A^)} 


( 11 ) 
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and that when lpa is normalized to be a probability density on [0,^4], 

— 5a{x)/A^ , where 5a{x) = 1{0 ^ x < A]{{x + 1) a (A — x + 1)]. (12) 


We will also make repeated use of a coupling of Komlos-Major-Tusnady [KMT76, Theorem 1]: 
let be i.i.d. random variables with zero mean, unit variance, and hnite exponential 

moments. For any A > 0 there are constants Ki,K 2 , and a coupling of to i.i.d. 

standard Gaussian random variables such that 


P( max 


- r.) 


2=1 


> iFl log fc + X < 7^26 


^—Xx 


3.1. Estimates for short time scales. We begin with some estimates for exp-minus-one 
walks run for k steps where k is arbitrary. We will apply these estimates for the case k < 
in the next subsection we derive better estimates for longer time scales k> A^. The following 
lemma is well known, see e.g., [PP95, Lemma 3.3]. 


Lemma 3.1. Consider a random walk whose step distribution has zero mean and unit 
variance. The probability that the walk started from x survives at least k steps before going 
negative is — {x + l)/{k + 1)^'^^, uniformly over k ^ 0 and 0 < x < 


Next we review an easy estimate for the probability the random walk stays conhned in an 
interval; this will be substantially rehned later for exp-minus-one walks. 


Lemma 3.2. Consider a random walk whose step distribution has zero mean and unit 
variance. The probability the walk will survive for at least k steps before exiting [0, A\ is 
< exp{— uniformly over all A > 1, k ^ 0, Xq e [OA]. 

Proof. Divide the time interval [0, k] into length-f subintervals with t ^ A^. It is then enough 
to note that the probability for the walk X to survive over a single subinterval is bounded 
away from one, uniformly over A and over the choice of the starting point x e [0, A\. Indeed, 
the survival probability is upper bounded by P(|Xt — Xq\ < ^4), which is bounded away from 
one either by direct calculation with the gamma distribution, or by applying the central limit 
theorem for (W ~ ^)/V^- Gl 

From now on we restrict our attention to exp-minus-one random walks in intervals [0, Al\ 
with absorbing boundaries. All our estimates hold also for one-minus-exp random walks (up 
to constant multiplicative error). Recall from (12) the dehnition of 5a- We assume from now 
on that A is at least of large constant size. 

Corollary 3.3. There is an absolute constant rj > 0 such that for k < {rjA)^, the probability 
for the exp-minus-one walk started from x to survive at least k steps before exiting [0, A\ is 
~ (5^(x)/(fc -t- 1)^/^, uniformly over 0 < 5 a { x ) < kfl‘^. 

Proof. As before it suffices to consider k exceeding any large fixed constant, so that k ^ k 1. 
Assume first that x < A/2, so Sa{x) = x -I- 1; the case x ^ A/2 follows in a symmetric fashion. 
From the one-sided bound of Lemma 3.1, the probability for the walk started from x to 
survive at least k steps before going negative is ~ (x -I- 1)/A:^/^, which trivially implies the 
upper bound. For the lower bound it suffices to subtract the probability that the walk 
exceeds A before going negative, which has probability < (x -I- 1)/A ^ p(x -f l)/k^/‘^. Taking 
sufficiently small p gives the lower bound. □ 
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Let p\{x, y) denote the probability density for an exp-minus-one walk X to go from x to y 
in k time steps without exiting [0, A\. Let (p\{x,y) denote the density at = y conditioned 
on survival in [0, A] for k steps. 

p\{x, y) = lim,|o {2e)~^¥{\Xk - |/| ^ e; X* e [0, A] for all 0 < t ^ A; | Xq = x), 
p\{x,y) = lim,^o (2e)"^P(jXfc - i/j ^ e | X* e [0,A] for all 0 < t ^ A;, Xq = x). 


Lemma 3.4. It holds uniformly over A, k, x, y that 


p\{.x,y) < 


5A{,x)5A{y) 

{k + 1)3/2 • 


Proof. As before it suffices to consider k exceeding any large fixed constant, so that A; — A; + 1. 
Take exp-minus-one walks X and Z started from Xq = x and Zq = 0 respectively, and take 
a one-minus-exp walk Y started from Yq = y, with X, Y, Z mutually independent. Fix also 
c e (0,1/2). The probability for an exp-minus-one walk to go from x to [y — e,y + e] without 
exiting [0, A] is upper bounded by the probability of the intersection of three events: 

El = {X survives at least ck steps without exiting [0, A]}, 

E 2 = {Y survives at least ck steps without exiting [0, A]}, 

FIs = {Z goes from 0 to Fcfc — ^ck + in (1 — ‘^c)k steps}. 

Now take c = 1/3, so the one-sided bound of Lemma 3.1 gives P(i?i) ^ 5^(x)/A:^/^ and similarly 
P(i?2) dAiy)/k^^‘^- The CLT then gives ^{E^lEi n — ^ck) ^ efk^A^ uniformly over 

all choices of T/fc ~ -Acfc- Multiplying these probabilities together proves the bound. □ 


Corollary 3.5. For k 

p\{x,y) 


A^, it holds uniformly over x,y that 

6A{x)6A{y) 


and PAi^^y) 


^Aiy) 


E, ^ 


A3 ^2 • 

Proof. Take the same notation as in the proof of Lemma 3.4, and take a constant c e (0,1/2) 
such that A^ < cA: < (i/A)^. Corollary 3.3 then gives P(Fli n E 2 ) “ 6A{x)6A{y)/k. Let 

_ j Z goes from 0 to T/fc — Xcfc + [—e, e] in (1 — 2c)k steps] 
without exiting the interval [—A/3,2A/3]. j 

Conditioned on Ei n E 2 , the event E^ that both Xck,Yck are in [A/3,2A/3] occurs with 
probability ~ 1: for 6a{x), dA^y) 3> log k this can be deduced directly from the KMT coupling; 
otherwise one can use the argument from Lemma 3.1 to reduce to the case 6a{x), dA{y) log k. 
The functional CLT gives 

P(E4 I El n ^2 n Eg) e/{k - 2ckYl^ ^ e/k^!^, 

so altogether p^(x, y) > 5^(x)5a(i/)/A:3/^ — 5a{x) 5A{y) / A^ . Combining with the upper bound 
of Lemma 3.4 proves PAi^^v) ~ (5a(x)5a(i/)/A3. A final application of Corollary 3.3 proves 
the estimate on the conditional density, ip\{x,y) ~ 5A{y)/A^. □ 


In preparation for the lemma that follows, observe that the Gam(t) density ft satisfies 


ft(t + x) ^ l/|x 4- 1| for alH ^ 1, X ^ —t. (13) 

For t = 1 it is easy to verify that /i(l -I- x) = ^ l/|x -L 1|, uniformly over all x ^ — 1. 

For t ^ 2, note that ft(t -T x) is unimodal in x, with the unique mode at x = —1. Since the 
density ft integrates to 1, for any x ^ —t we have 

rh+i| 

ft(t - 1 + u)du ^ \x + l\ft(t + x), 

Jo 
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Lemma 3.6. It holds uniformly over A,k,x,y that 

(rr < J ^ + i)/l^ ~ y\^ ^ 

Pa\ ’ ~ I _ X + 1)(|/ -f l)/jx — for X > y. 

Proof. Assume 0 ^ x < y ^ A\ the case x > y follows by a symmetric argument. Define 
x' = X + {y — x)/i, y' = y — {y — x)/4; note that we need only consider y — x exceeding a large 
constant. Take exp-minus-one walks X and Z started from Xo = x and Zq = 0 respectively, 
and take a one-minus-exp walk V started from Yq = y, with X, Y, Z mutually independent. 
Define the stopping times a = minjt ^ 0 : ^ [0, x'f] and r = min{t ^ 0 : Yt f [y', A]}. The 

probability for an exp-minus-one random walk to go from x to [y — e,y + e] in k steps without 
exiting [0, A] can be upper bounded by the probability of the intersection of three events, 

E, = {X, > a:'} , 

E, = {Y,<y'}, 

E^ = {a + T ^ k} n {Zfc_(a-+T) ^ ~ f]} • 

It is clear from Lemma 3.2 that the stopping time a is integrable, so Wald’s identity gives 
X = EXo-. Rearranging gives 


P(T^i) = --where a = EiX^. I X„ < 0] and b = E[Xo-1 X„ > x'l. 

b — a 

If X is either the exp-minus-one or one-minus-exp walk then — 1 ^ a ^ 0 and x' ^ b ^ x' + 1, 
so P(i?i) ^ (x -I- l)/x' < (x -I- l)/{y — x), and similarly P(i? 2 ) — V + ^)/{y — x)- Condi¬ 

tioned on El (A E 2 , we have y' — 1 ^Yt- ^ y' a.s., while X„ — x' is a standard exponential 
random variable. For X„ — x' large we take the crude bound 


F{E 3 -,Xcr — X > [y — x)/4) ^ max max 

t^k ugM. 


F{Zt e M + [—e, e]) 


^{y-x)/4 


< 


< 


g(y ~ y 


X 


For smaller values of X„ — x' we instead bound 


P(Fl3; Xa- — x ^ (y — x)/4) ^ max max F(Zt e (y' — x) — u + [—e, el). 

t^k 0!S:U^l+(y—x)/4 

which is again < e/{y — x) by (13). Multiplying the probabilities together gives the claimed 
bound. □ 


3.2. Estimates for long time scales. Recall from (11) and (12) the definitions of Aa, Pa-, Sa- 
Lemma 3.7. For k > A^ , it holds uniformly over x^y that 

PAix,y) ^ and (p\{x,y) ^ 

Proof. Take t ^ k with A^ < t ^ (hA)^, and run the walk for t time steps: by Corollary 3.3, the 
walk survives in [0, A] up to time t with probability “ 5a{x)/A^'^ ~ 5a{x)/A. By Corollary 3.5, 
the density at X^ = 2 ; conditioned on survival is (p\{x, z) ~ 5a{z)/A^. In particular we can 
hnd a large absolute constant C such that 

C'"Va(^) < ^ C^a{z). 

Therefore, if we start from the distribution (^a(x, •) and evolve the walk forward for A: — t steps, 
the terminal density at X^ = y (conditioned on survival up to time t) will be sandwiched 
between C'“^(Aa)^“Va(i/) and C(Aa)^“Va( 2/)- Both lower and upper bounds agree up to 
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constant factors with {\aY 5Aiy)/, and multiplying with the probability of survival up to 
time t proves the hrst estimate y) “ {XA)’^SA{x)6A{y)/A^. Integrating over y proves that 
the probability of survival up to time A; is — {Xa)^Sa{x)/A, therefore y) — 6A{y)/A'^. □ 

Remark 3.8. Since p\{Aj2, Aj2) is nondecreasing in A, an immediate consequence of 
Lemma 3.7 is that is nondecreasing in A. 


Recall that (A^)^ 

p\{x,y) < 


1 for k < A^. It therefore holds uniformly over all A, k, x, y that 
(Ayi)^^A(x)5yi(|/) by Lemma 3.4 for k < A^, 


(A^ A (fc + l))3/2 

For A — X — y ^ A we also have the bounds 

p’Xix^A-y)] ^ (Aa)^(x +l)(i/ +1) 


and Lemma 3.7 for k > A^. 


(14) 


p\{A-x,y) 


< 


A3 


by Lemma 3.6 for k < 
and Lemma 3.7 for k>A^. 


(15) 


In both (14) and (15) the < can be replaced with — in the regime k> A^. 

We conclude with our estimates for the exp-minus-one fc-bridge, which we recall is an 
exp-minus-one walk conditioned to return to the origin at time k. We hrst estimate the 
probability R\ that this process has range at most A, 

P(| range(X)| < A; \Xk\ ^ e) 


Ra = lim 

elO 


P(|Wfc| < e) 


X an exp-minus-one walk. 


Lemma 3.9. For k > A^, the exp-minus-one k-bridge has range at most A with probability 

^^3/2 


Ra - M 


A3 


Proof. Let X be an exp-minus-one walk started from Xq = 0. Decompose R\ = R\i^ + 
where Ra /4 contribution from the event that the range of X is smaller than A/4, while 

R}a ^ is the contribution from the event that the range is in [A/4, A]. In order to have range 
less than A/4 the walk must certainly stay conhned within distance A/4 of the origin, so 


Ra/a < lim 


^A/i 


elO 


P((AWs|^(|V4V4];|.U.|^e) ^ < k^'\\Ant1A 


by (14). For k > A^ the right-hand side is < k^l‘^{XAY/A^. 

For larger H we make a more precise calculation. Let S,T e [0,A:) denote the times where 
the walk achieves its minimum and maximum respectively: the contribution to R\ > from 
5 < T is 


^ P(range(X |[o,,]) EdH;S< T; |Xfc| 


[ 

Ja/4 

k-l 


^ e 


P(|X,| < e) 


dH 


k—1 ^3A/4 

k^R 2 l{s < t} p)7^(0, H)p%-^+^{H, 0) dH < k^R 
s,t=o 4 o 


(Aa- 


(A-A)' 


dX 


where in the hnal step we applied (15) and made the change of variables H = A — A. From 
the expansion (11) and the monotonicity of Aa (Remark 3.8) we have 


A 


A-A 

"aT 


^ max 11, exp 


-0(A)-L 0(1) -I 
A3 J 


(16) 











SUPERCRITICAL MINIMUM MEAN-WEIGHT CYCLES 


11 


SO we find that the contribution to R\ ^ 


<k^l^ 




J 

Jo 


from S' < T is 

3A/4 


^A-A 


dA - 


(^aY 

A3 


The contribution to ^ from S > T has the same value, so the upper bound follows. 
The lower bound can be obtained in a similar manner, but summing only over pairs s < t 
with t — s > and k — t + s > A^ and applying the lower bound from Lemma 3.7 (see the 
comment below (15)). □ 


We conclude with an estimate on the local time profile for the exp-minus-one bridge, which 
will be used to control the variance of “good” cycles in our second moment argument. Let X 
be an exp-minus-one fc-bridge, and let X be X with a constant shift that centers it at A/2\ 

Xt = Xt +{A- maxX - minX)/2 (17) 


We then dehne the local times of the range-restricted /c-bridge by 


tAS) ^ E 


^llWeS'} |range(X)|^A 


i=l 


for S c [0,T]. 


Lemma 3.10. For k > A^, it holds uniformly over subintervals S c [0, A] that the loeal time 
of a range-restricted k-bridge is bounded by 

(YYS) < (1 + \S\) maxdAixy. 


Proof. Let X denote an exp-minus-one walk, and H its exact range, so that the recentered 
walk X is conhned in the interval Ia,h = [(^ — + -^)/2]- We will compute 


^\iS) = lim 

ej,0 


j:l,F{Xte S,H ^A;\Xk\^e) 

F{H ^A-,\Xk\^e) 

e{\Afk/A3 (since X is a walk rather than a bridge. 


Lemma 3.9 gives that the denominator is 
there is a factor of e/\/fc)- We decompose = separating the numerator into 

the cases H < A/A and A/A ^ H ^ A. For H < A/4, applying the bound (14) gives 


4<(^) 




4<([o,^]) 


< 


^Pa/2(^/4,^/4) 

{^AYk/A3 


^^{^A/2/^aY ~ ^ 


< max 5 a (t)' 

xeS 


(18) 


where in the last bound we used that X is confined in pA/S, SA/S], so in fact is zero 

unless 5a{x) — A for some x e S. 

For larger H, we sum over all possibilities s,t e [1, /c] where the walk achieves its minimum 
and maximum respectively, as well as all times u e [1, fc] where X„ e S. For simplicity we 
consider only the contribution from times s < u < t] the contribution from other permutations 
is calculated similarly and will be of the same asymptotic order (up to constants): thus 

d»(S)S L // A/A r f p'‘„-‘{0,x„)p‘i;'‘{x„,H)p%-‘*‘{H,0)dxdH 

s,u,t=l JA/4JsnlA,H 


where xh = x — [A — H)/2, so that for x e Ia,h we have xh ^ [0, H] and Sh{xh) < 5a (t). By 
symmetry it suffices to consider xh ^ H/2: then pff^{0,XH) can be bounded by (14), while 
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H) and can be bounded by (15). The total contribution to 

from times s < u <t with u — s ^ is 

while the total contribution from times s < u <t with u — s ^ is 

Adding these together and combining with (18) gives the stated bound. □ 

4. Light cycles and long cycles 

Recall that for a given cycle ^ we use c('^) = n wgt(‘^)/len(P) to denote its mean weight 
scaled by n. In this section we prove that w.h.p. no cycles in the supercritical regime c> Ije 
have Co — c 3> l/(logn)^. We also show that in the regime |c — Co| l/(logn)^, w.h.p. there 
are no cycles of length (logn)^. The formal statement is as follows: 

Theorem 4.1. In the complete graph or complete digraph with i.i.d. unit-rate exponential 
edge weights, for all e > 0 there exists a constant C = C{e) > 0 such that 

P (3 cycle "to : 1/e < c{f€) ^ Co — C/(logn)^) < e. (19) 

For all e > 0 and all Ci > 0 there exists a constant C 2 = C 2 {e, Ci) > 0 such that 

3 cycle : 1/e < c(‘^) ^ Co + Ci/(logn)^ 
and len(‘^) ^ C' 2 (logn)^ 

4.1. First moment for uniform cycles and paths. We begin with some hrst moment 
estimates for uniform cycles and paths. Let 

Z^{A) = number of c-light A-uniform fc-cycles; 

Z^{A) = number of fc-paths that, for some A" e [A — 2, A], (21) 

have c-excedance —A" and are (c, A")-uniform. 

We will also denote 

ZpW - L Zi(A), Z,(A) = Z>‘^(A), ZJA) ^ 21 A'O- 

The purpose of defining Zc(A) is to bound the number of light cycles which do not contribute 
to Zc(A). From Dehnition 2.3, if a cycle ^ is c-light but fails to be A-uniform, then it has 
a subpath P with c-excedance < —A. The following lemma shows how to extract further 
subpaths from P that are both light and uniform: 

Lemma 4.2. For 2 < A' < A, given a path P with c-excedance < —A, one can extract at 
least 1 -L [A — A'J distinct (though not necessarily disjoint) subpaths p ^ P, each contributing 
to Zc(A') as defined in (21). 

Proof. By scaling we can assume c = 1. Suppose P has length k and edge weights Wi,... ,Wk. 
By assumption, the process Wj = '^l^finwi — 1) goes from ITo = 0 to Wk < —A. Let 

{xi, yt) = -{i,i + A'), for 0 < i ^ [A - A'J. 

Since A' ^ 2, p, -f 1 < a:* — 1 for each i. Since W decreases by at most one at each step, for each 
i there is some pi for which 11/,- e (xj — 1, Xj], and some Tj > p, for which R/-. e (p*, pj -I- Ij. Let 
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\fji, fi\ be any minimal subinterval of [ 0 , k\ for which Wfn e (xj — 1 , Xj] and Wf^ e (?/j, + 1 ], 

and let Pi be the subpath of P corresponding to the interval \fii,fi\. Since is minimal, 

the subpath pi is {Wf^^ — bhfj-uniform, where A' 




X 


Wf < A'. Since the intervals 

‘ i ^ 

□ 


l,Xj] are disjoint, the paths Pi are distinct. 

Applying Lemma 4.2 with A' = A, when A ^ 2 , we see that any c-light cycle either 
contributes to Zc(A), or has a subpath contributing to Zc(A). 

Lemma 4.3. For k > A^ , c < 1, and 0 < 5 < 1, 

{n)k (ceA,i)‘[l - (1 - i)‘] 


EIZ;(A) - Zy_„(Al)] 


( 22 ) 


nk ^3 

Proof. Recalling Dehnition 2.3, a fc-cycle with edge weights wi,... ,Wk summing to kc/n 
contributes to Z^{A) if and only if c < c, and the untilted bridge with increments {nwifc — 1 ) 
has range at most A. Applying (7) with s = c, we see that EZ^(A) = (EZf)R\, and then 
(22) follows from Lemmas 2.1 and 3.9. □ 


Lemma 4.4. For A ^ 2, c < 1, and k ^ 1, 


EZ^{A) < 


{ce\AY ^ 


A3 e 


A' 


(23) 


Proof. A fc-path with edge weights wi,... ,Wk summing to kc/n contributes to Z^{A) if and 
only if, for some A" e [A — 2 , A], the process with increments {nwi/c — 1) goes from 0 to 
—A" without exiting [—A",0]. In particular this implies c = c(l — A"/k). Making the change 
of variables A" = A — A, we have 


I 


EZ^.iA) = in)k+i fk[{k - A + A)c/n] 


Pa_^{A — A,0) cdA 


to fkik - A + A) n 

where the hrst factor in the integrand is the probability density of the path weight, and the 
second factor is the probability that an exp-minus-walk, conditioned to go from A — A to 0 
in k steps, does so without exiting [0, A — A]. Simplifying gives 


n' 


I 


EZ;(2l) « I - A,0)dA 


< 


n [ccXaY 
^ A3 


Aa-A 


dA, 


where the last bound follows from (15). The final integrand is < 1, which gives (23). □ 


4.2. Uniform cycles and light cycles. We now apply our hrst moment estimates to argue 
that in the supercritical regime c > 1/e, cycles cannot be too uniform or too light. Recall 
(10) that we set Aq = logu and Co = l/(eAAo)- More generally, we dehne the relation 

ccAa = expj-y/A^} (24) 

— where we can think of c as being determined by A and 7 . A hrst consequence of Lemma 4.3 
is that cycles of length at least (logu)^ cannot be too uniform in the targeted regime for the 
mean weight: 


Lemma 4.5. For all e,C > 0 there exists A = A(e, U) > 0 large enough that 


P 


3 cycle ^ : c{f^) < Co + U/(logn)^, 
len(^) ^ A^, and ^ is A-uniform 


< e 


for A = logn — A. 
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Proof. Suppose (24) holds with 7 > 0. Using Markov’s inequality with Lemma 4.3 and (24), 
F{Zf^\A) > 0) < EZf^\A) < < I /7 , 

k^A'2 

SO by taking 7 ~ 1/e we can ensure P{Zf^^ [A) > 0) ^ e. Using the expansion (11) for A^i, 
we see that by choosing A = logn — A with A — U + 7 , where 7 and C are hxed as n —>■ 00, 
we can ensure c = l/(eAA) exp(— 7 /A^) ^ Co + (^/(logn)^. □ 

We next show how Lemma 4.3 and Lemma 4.4 can be combined to rule out cycles that are 
lighter than the targeted regime for the mean weight. 

Proposition 4.6. For all e > 0 there exists C = C{e) > 0 such that 

™ / 3 cycle : ciho) ^ Co — (^/(logu)^ 
y and len('^) ^ (logu)^ 

Proof. From Lemma 4.2, assuming A ^ 2, a c-light fc-cycle either contribntes to Z^{A), or 
has a snbpath contribnting to Z^{A) for some i ^ k. By Lemmas 4.3 and 4.4, if we take 
A = Ao and 7 a large positive constant, then 

p--ye/A^ 

P(Zf > 0) < + HZ,(A) S L = l/T . 

so by taking 7 ~ 1/e we can ensnre > 0) ^ e. Since c/co = exp(— 7 /yl^) and Co ~ 1, 

the proposition follows with U ~ 1 /e. □ 



4.3. Non-uniform cycles and long cycles. Recall that in Lemma 4.5 above we ruled out 
cycles that are too uniform in the targeted regime for the mean weight. We now prove the 
complementary assertion that cycles cannot be too non-uniform in this regime: 


Lemma 4.7. For all e, C* > 0 there exists A = A(e, U) > 0 large enough that 


P 


3 cycle : c(^) ^ Co 3- CIffogn'f \ 


^ e. 


and ^ is not (logn -t- A)-uniform J 

Proof. We can set 7 = 1 and A = log n — A' with A! C and have c = l/(eA^/) exp(— 7 /(A')^) 

Co 3- U/(logn)^. Then Lemma 4.4 gives 

n 1 


E^c(A') < 


n/e^ = 


{A'fe^' 1 - e-AA')^ 

Let A = logn 3- A ^ A'. By Dehnition 2.3 and Lemma 4.2 together with Markov’s inequality. 


P(Zo > Z,(A)) < P(Ze(A') ^ 1 + [A - A'J) 




EZo(A') 


< 


1 + [A 3- A'J ~ (A 3- A') 

By taking A ~ e^'/e we can ensure that this probability is sufficiently small. 


□ 


In the remainder of this section we argne that there are no cycles in the targeted weight 
regime that have length much longer than (logn)^. In view of Lemmas 4.5 and 4.7, it remains 
to consider cycles which are neither too nniform nor too non-nniform: 


Lemma 4.8. For all e,C, A > 0 there exists B = B{e,C,A) > 0 large enough that 


/ 3 cycle ^ : Co — C/(logn)^ ^ c(^) ^ Co 3- C/(logn)^, 
F \ ^ is (logn 3- A)-uniform but not (logn — A)-uniform, 
y and len(‘^) ^ i?(logn)^ 




e. 
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Proof. We denote 

c-= Co + C / ifognY, A- = logn + A. 

Consider a fc-cycle with weights Wi,... ,Wk summing to kc/n, where c“ ^ c < c"^, which is 
A’^-uniform but not y4“-uniform. Let X denote the corresponding untilted bridge, with 
increments {nwijc— 1). For convenience we shift X by a constant so that its minimum is 
zero. Then its maximum is in [y4“,y4''']. 

Let H denote the ceiling of the maximum of X. Write H = [Ao] + Xh (so Xh is integer¬ 
valued). Let A' = [Ao] — X' be the largest integer for which A' < A~ — 2 and 

c'^cXa < exp{ —l/A^} for all A such that (logn)/2 ^ A ^ W -I- 2. (25) 

The constraint A' < A~ — 2 ensures A' -L Xh ^ 1. It follows from (8) that A' ~ max(C, A). 
Define 

X = H - {A' -2) = Xh + X' + 2, 
y = A - 2 = H - {Xh + X' + 2). 

For fixed C and A when n is big enough we have 0 < x < y < H. 

For a sequence of times 0 < ti < ... < tg ^ k with q even, let E\ti.,g\ denote the event that 
the ti partition the bridge X into alternating up-crossings and down-crossings: that is. 


E\ti.,g\ = {U\ti_iAi\ for all i even} n {D[L_i,L] for all i odd} 

where [/[a, h\ and D[a, h\ indicate the events of up-crossing and down-crossing on [a, h\\ 

f/[a, h\ = {Xa < X and b = min{f > a : Xt > y}} 

D[a, b] = {Xa > y and b = min{f > a : Xt < a:}}. 


(with cyclic time indexing). The proof idea is as follows: any cycle with a very large number 
of crossings will give an excessive contribution to Zc{A") where A" = y — x, so we can assume 
a bounded number of crossings. We then argue that a long cycle with a bounded number of 
crossings is unlikely to occur because its range is effectively reduced from A to A. 

Each down-crossing corresponds to a subpath with c-excedance < —A". By Lemma 4.2, 
from each such down-crossing we can extract a contribution to Za+{A"). From Lemma 4.4 
and (25), the expected number of cycles with at least Q/2 down-crossings is 




¥.Zo+{A‘ 

Q/2 


< 


n 


Qe 


A" 


L 


(c+eAyi")^ 


{A") 


< 


n 


,,A+2A' 


Qe 


A" 


< 


Q 


which can be made arbitrarily small by taking Q large. 

It therefore remains to control the cycles with q < Q. To this end, we shall first bound the 
probability for an exp-minus-one walk Y started at Yq = Xq — with no conditioning on the 
value of Yk — to belong to the event E{ti.q\. 

1. First consider the probability of a length-f down-crossing, D[a,a-ff], conditioned on 
having just completed an up-crossing at time a. This means that [Wl = y + h for some 
integer 1 ^ h ^ H — y. Conditional on h, the probability to make a down-crossing is 
upper bounded by the probability that an exp-minus-one walk started exactly a,t y + h 
will travel to [x — 1, a; -t-1] in t steps without exiting the interval [x — 1, if -t- 1]. It follows 
that the probability of down-crossing D[a, a + t], conditioned on having just completed 
an up-crossing at time a, is upper bounded by Pdn{t) where 

Pdn{t) = max p\,{y + h — {x — l),u) du. 

IsihsiH-y Jg 
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Recalling H — y = Ah + A' + 2 and applying (15) gives 


Pdn{t) 


< 


{Xa'Y{Ah + A') 


2. Similarly, consider the probability of a length-t up-crossing U[a,a + t], conditioned on 
having just completed a down-crossing at time a. This is upper bounded by the probability 
Pup(t) that an exp-minus-one walk started exactly at the integer x will travel to [y, H + 1] 
in t steps, and remain in the interval [0, i/ -t- 1] c [0, A'] up to time t — 1. Conditioning on 
the position y + 1 — u of the walk at time t — 1 and the position ^ at time t gives 


Pnp{t) 


< 


rH+l r 

Jv Jo 


p\Y{x,y + 1 - u) 




dudz < 


p\Y{x,y + I - u) 


du. 


If t > (logn)^, then (14) gives 


^ {Xa'Yx r+i u + l {Xa')\^h + A') 

~ (A03 Jo e“ (A')3 

If t < (logn)^, the contribution from u < y/2 satishes the same bound by applying (15) in 
place of (14), while the contribution from u ^ y/2 is negligible in comparison. 



Conditioned on all these crossings, the probability that jlfc — ho| < e is < e, whereas uncon¬ 
ditionally we have \Yk\ < e with probability “ e/k}/"^. Altogether we hnd 

F{X e Eli.J) - lim-p(|y^| ^ ^-~ -(1^^-' 

Therefore, the probability to have a cycle with |c(‘^) — Co| ^ C/(logn)^, which is A"*"- 
uniformly light but not A”-uniformly light, has length len(^) ^ L = i?(logn)^, and q ^ Q 
is, using Markov’s inequality and Lemma 2.1, 


L PA E £[«,„])< 2 ] 


{c+eXA')^ sp + A')A;]''eO('?) 


k^L 


k^L 


km 


T 




We can assume B > 1/{Ah + A') to make the summand increasing in q. Recalling (25), this 
in turn is 


< 


(A')3(Aj, + ^ Q[k/{A'f]^ 


B{\ognY 


k-^L 


exp{fc/(A')3} {A'Y 


< 


Q{Ah + r 


B 


JB 


— dz. 


This can be made arbitrarily small by taking B large, and the result follows. 


□ 


Proposition 4.9. For all e,C > 0 there exists B = B{e,C) > 0 large enough that 

™ / 3 cycle ^ Co — C/{\ognY < c('^) ^ Co 3- C/{\ognY 
y and len('^) ^ R(logn)^ 

Proof. Follows by combining Lemmas 4.5, 4.7, and 4.8. □ 

Proof of Theorem f.l. Proposition 4.6 rules out light supercritical cycles of length ^ (logn)^, 
and Lemma 2.2 rules out the ones shorter than (logu)^, which implies (19). Proposition 4.9 
rules out cycles of length 3> (logn)^ in the regime |c — Co| ^ l/(logn)^, which together with 
(19) implies (20). □ 
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5. Variance of typical-profile cycles 


In this section we complete the proof of our main result Theorem 1.1, the key ingredient of 
which is to demonstrate that cycles of length ~ (log n)^ exist in the regime |c — Co | ^ l/(log n)^. 
(We also argue that it is unlikely that much shorter cycles exist in this weight range.) We 
prove the existence of these cycles by a second-moment computation on a set of “good” cycles, 
as we now formally define. 

5.1. Good cycles. Recall Ao = logn, and let A = [Ro] — A for A a large positive integer, 
with A ^ (logn)/2. Choose c = ca such that ccXa = 1- 

Definition 5.1. A fc-cycle with weights Wi, ... ,Wk summing to kc/n, is termed A-good if 

(i) c(l — 1/k) ^ c ^ c ioT c = Ca, A = [Aq] — A; 

(ii) the process X with increments {nwijc — 1) has range < A — 2, 

(hi) the recentered process X (as in (17)) has < A x ^^(a;)"^ visits to {x — l,x], for each 
1 < X < A. 


Let flA{k) be the number of A-good fc-cycles and let 

Da = ^ DA(fc). 

fc=A3/A 

We show by the moment method that for large A, Da is positive with large probability. To 
this end we first argue that Da is large in expectation: 


Lemma 5.2. For A large, ED a ~ A. 


Proof. For A^ < k < consider a fc-cycle ^ that satishes properties (i) and (ii). The 
conditional probability that (iii) fails is, by taking a union bound over 1 < x ^ A and applying 
Lemma 3.10 with Markov’s inequality. 


< 


L 

X=1 


A6A{xy 


< 1/A. 


It follows, assuming A large enough, that the expectation of Da(/c) is equal up to constants 
to the expected number of fc-cycles satisfying properties (i) and (ii) only. By Lemma 4.3, 

ESlAk) » E[Z^(A) - Zy_v«(-4)] = 7^ _ A. 

Summing over A“^ ^ k/A^ < A proves the claim. □ 


5.2. Variance bound. The crux of the proof is in the following lemma which argues that 
the count of good cycles has low variance. 


Lemma 5.3. For any cycle 

P(^' is A-good\‘T^ is A-good) ^ EDa + oa(1) 

where oa(1) indicates an error tending to zero in the limit o/ A — > oo. 
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Proof. Fix A large, and fix a k-cjcle Condition on the event that ^ is good (from here 
on, “good” means A-good), and write 

S(^) = ^ P(^' is good I ^ is good). 

If a cycle is disjoint from ^ then the events is good} and is good} are independent: 
thus, the contribution to from all n = 0 is at most EIIa. It remains then to 

argue that the contribution from cycles intersecting with ^ is oa(1). 

To this end, consider ^ ^ such that ^ n consists of q shared segments, with g ^ 1. 
Let r count the vertices in . The cycle is partitioned into alternating in-segments 
(contained in and out-segments (edge-disjoint from ^); the combined length of all the 
out-segments is r -I- g. Label the in-segments \uj,Vj\, 1 < j ^ g, in order of their traversal by 
^ may traverse the segments in a different order. Let nout(L, g) count the number of 
such cycles for a given tuple (f^,r,q, {uj,Vj)i<ij<iq), 

nout(r, ^ ^ I X ^ ^ 

(since is given by choosing an ordered sequence of r vertices, then dividing them in to g 
groups). 

Let X' be the process with increments {nw'Jc — 1), where w\ are the edge weights on 
The edge weights in are still i.i.d. unit-rate exponential random variables after the 

conditioning on ^. For 1 ^ j < g let integers yj, Zj be defined by 

yj = (discretized increment over j-th in-segment); 

Zj = (discretized increment over j-th out-segment). 

Write y = {yi,... ,yq) and likewise z = {zi,..., Zq). Let y denote the average of the yj, and 
likewise z the average of the Zj. For to satisfy property (i), we must have q(fy -\-P) e { — 1,0}. 
For each j, yj A ^ 1 and —Zj -I- A ^ 1. Thus, given y, the number of choices for compatible 
(y, z) is at most 

Next, by property (iii), the number of choices of {uj,Vj) consistent with y is at most 
n{A^J]5A(a:)^5A(a:+ |%|)^} ^ A^*? |(A - |i/^|) max [(x + 1)7^ - |i/j | - a; + 1)^] |, 

i=i X j=i 

where we have used that there are < A — \yj\ choices of x for which the summand is positive. 
For each j the maximum is achieved with a; -I- 1 = (A — \yj\)/2, and combining with Jensen’s 
inequality the above is 

< f[(-4 - I^D’ « A"’(^ - V' L l!/,l) < + y)'» = b,®. 

i=i j 

Lastly, let p(z) be the probability that X' has increments as specified by z, with range at 
most A (otherwise 'if' would violate property (ii)). Writing ij for the length of the j-th 
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out-segment, a similar calculation as in the proof of Lemma 4.4 gives 

^ I I f (J{cjn) f ( ^ ^— ma.xpj^[x,x + Zj) 




fi Ai + Zj) 


Applying (15) and making some straightforward manipulations (recalling ccXa = 1), we arrive 
at the conclusion that for any 2 consistent with y, 


P{z) ^ 


ceX, 


n 


c^^^AcXa (A — z) 


^ 12 -. 


ne~ 


A3 


^ e 


0{q), 


^r,q V 


Cr,qiy) — „ 


1 r(A + ^)2^ 


rf L nevA^ - 


Combining these estimates gives 
S(^) - Ef^A < ^ nout(r, g)ag(y)bq(^)c 


r,q V 






rgO(i)^ 2 (A + yY^-tq{r + q) 


9-1 


nev 


A39 


qx,y qx,y 

where the sum is taken over g<r + g< AA3, and \y\ ^ A — 2 with qy integer-valued. 

Making the change of variables Y = A -I- y ^ 2, we hnd 

q 


S(^) - EQa ^ Yj [ 


r A3 yi 3 . 

9 ^ rA3yi3 

Vn/e^ eX . 

q,Y 


q,Y 

which tends to zero in the limit A —> oo as claimed. 
Corollary 5.4. In the limit A —> oo, P(f2A = 0) —> 0. 
Proof. By Chebychev’s inequality and Lemma 5.3, 


A3 


< 


□ 


p(riA = 0) ^ 


Var|f2A| EIL^aI + oa(1) 


< 


(E|f2A|)2 ■ 

By Lem. 5.2 the right-hand side tends to zero in the limit A —> oo. □ 

5.3. Conclusion. We conclude this section with the proof of our main result. 

Proof of Theorem 1.1. Theorem 4.1 and Corollary 5.4 together imply that we can choose 
Cl = Ci(e) sufficiently large so that 

Co - Ci/(logu)3 ^ nWn < Co -t- C'i/(logn)^ 
and (logn)^/Ci < < C'i(logn)3 

Combining this with Lemma 4.3 and Lemma 4.7, we can choose C 2 = C 2 A, Ci) sufficiently 
large so that for any interval I c [ 0 , Ci] with |/| ^ l/C' 2 , 

3 cycle : [c — Co](log?T,)3 e I 
or len(‘^)/(logu)3 e I 

Combining these proves (2) and (3). □ 


P 


nWn > 1/e I ^ 1 — e. 


P 


^ e. 


6. Eigenfunctions of the one-minus-exp walk 

Recall that the one-minus-exp random walk is the real-valued stochastic process (W)t>o 
whose jumps Xt+i — Xf are independent and identically distributed as 1 — Exp(l). In this 
section we study eigenfunctions of this walk with killing outside an interval [0, H]. 
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6.1. Eigenfunctions and eigenvalues. Consider the one-minus-exp walk killed outside 
[0, i/], and suppose that is a left eigenfunction of this process with eigenvalue A: this means 
(see Remark 6.1) that is a function supported on [0,iir], not identically zero, satisfying 

rH 

Xg{x) = ^ l{u ^ X — l}g{u)e ^ du, 0 ^ x ^ H. (26) 

Jo 

As g is not identically zero, {x : g{x) 7 ^ 0} is a non-empty subset of [0, H], and we denote its 
inhmum by Xmin- We then have 

fH 

Xg{x)e^ ^ = g{u)e ^ du, 0 ^ x ^ min{i7, 1 -t- X min }. 

Jo 

— including for some Xmin < x < mm{H, 1 -f Xmin} for which g{x) ¥= 0. If the integral on the 
right-hand side vanishes, then A = 0; otherwise A 7 ^ 0, Xmin = 0, and (by rescaling) g{x) = 
for 0 < X < min{if, 1}. In particular, for if < 1 we conclude the process has a unique non-zero 
eigenvalue A = Hje with associated left eigenvector g{x) = e^l{0 < x ^ H]. 

We see from (26) that g is smooth on (0, H)\L, with continuous derivatives up to order k 
on (0,ii)\Z^fc. We can therefore differentiate (26) to find 

Xg'{x) = Xg{x) — g{x — 1), 1 < x < ii. (27) 


Take for the moment H = cc. We can solve for g on the intervals [k, A; - 1 - 1] {k e Z^i) one at 
a time, as follows: if g satisfies (27), then h satisfies h'{x) = g{x — 1)/A for all x > 1, where 


h{x) = —g{x) + J g{u) du, x ^ 1 . (28) 

Knowing g{x) for 1 ^ x < x determines h{x) for 1 ^ x ^ x; the reverse also holds since 

r 

0 = —g{x) + g{u) du, 0 < x ^ x 


has only the trivial solution g = 0. Suppose inductively that g (hence h) has been determined 
on [0, k] for k e then we can determine h on the next interval [A:, A: -I- 1] by evaluating 

rx—\ 

h{x) = h{k) - 1 - A“^ g{u) du, A: ^ x ^ A: + 1, 

Jfc-i 

which determines g on [A:, A: + 1]. This proves that in the case H = ao, there is a unique 
continuous function g = g\ supported on [ 0 , 00) that satishes g{x) = for 0 < x ^ 1 , and 
satisfies (27) for x > 1. It is straightforward to verify that 


9\{.x) 


Pi-i 


k=0 


(x — kY 

{—XeYkl ’ 


X ^ 0 


(29) 


is such a function, so it must be the unique one. 

If for hnite H and A 7 ^ 0, (7 is a solution of (26), then g satisfies (27), and it follows from 
the preceding discussion that g must be the restriction of gx to [0, 77]. We calculate 

1 . _u , _ (a -1- 1 — kY ^ (6 + 1 — kY _ g\{a + 1) gxY + 1) 

Xejjx[u)e (-Ae)^7! (-Ae)^A;! “ e“+i ’ 
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Figure 1. The left plot shows gi (x), and the right plot shows 50.95 (x) ■ 



Figure 2. The eigenvalues of the walk with two-sided killing as a function of H. 


SO equation (26) with g = g\ and 1 ^ x ^ H simplifies to 


0 = 


9\{x) 


1 

Ae 


r 

Jx—l 


' du = 


gx{H + i) 




therefore g\{H + 1) = 0 which is a polynomial equation in A of degree \H] (Fig. 2). 


Remark 6.1. We comment briefly on left and right eigenfunctions of the process. Formally, 
the operator for the one-minus-exp walk killed outside [0, H] is given by 

{Kf){x) ^ E.[/(Xi)l{0 < < H}] = f f{y)p{x,y) dy, 

Jo 

where p is the transition kernel from x to y. The left action of K is the right action of the 
adjoint operator, {K*f){y) = {fK){y). In the current setting, p{x,y) = k{y — x) where k is 
the density function of the 1 — Exp(l) random variable. We therefore have 

nH nH 

iK*f){x) = f{y)k{x -y)dy=\ f{H - y)k{y - (H - x)) dy = {RHKRHf)ix), 

Jo Jo 

where Rnfix) = f{H — x). The reflection Rh is involutive, and it relates the left and right 
eigenfunctions of 5 ^ is a left eigenfunction of K (i.e. a right eigenfunction of K*) if and only 
if Rng is a right eigenfunction of K with the same eigenvalue. If /i, /2 are right eigenfunctions 
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of K with non-zero eigenvalues Ai, A 2 then, writing (•, •) for the -h^[0, if] inner product, 
= iX2)-\RHfl,Kf2) = iX2)-\K*RHhJ,) = (X^/X^KRh 
Consequently, if Ai ¥= X 2 then (Rnfi, / 2 ) = 0. 


6.2. Series expansion. We now review the solution obtained by Wright [Wri49] for general 
homogeneous difference-differential equations with constant coefficients, 

m n 

A{y){x) = 2 Zj = 0’ e C, 0 = 60 < < ■ ■ ■ < (30) 

/i=0 u=0 

where denotes the v-ih. derivative of y. Observe that 




A(g^x) _ = Z Z 

/X = 0 !^ = 0 

so if r(s) = 0 then y{x) = solves A{y) = 0 . More generally, we can apply the formula 
= 2^=0 to calculate 

A(x^en = Z Z Z 

^=0 i^=0 r'^Q 

where for r e [j)r denotes the falling factorial j{j — 1 ) ■ ■ ■ (j — r -I- 1 ), which is zero if 
r — j E Expanding in powers of x gives 

e 


i=o 






fi=0 u=0 r^O 


T ! 


If s is a root of r of order £ -I- 1, meaning 0 = r(s) = t'{s) = ... = then A(xR^^) = 0 

for the integers 0 < j < 

Assume that m,n ^ 1, and that the m x n coefficient matrix contains a non-zero entry 
in each of the first and last rows and columns, which eliminates pure difference equations and 
pure differential equations. Under these assumptions, it is shown [Wri49] that the general 
solution to (30) is given by a limit of linear combinations of the solutions xR^^ described 
above; further, it is explained how to compute the coefficients of this linear combination 
given initial data 2 /*''^^( 0 ) (0 < z/ < n) and y^'^\x) (0 < x < bm)- Let 

m n 

Hi{s) = Z Z y^''\u)e-^'^du, 

At=l u=0 

m n u 

H(s) = Hi(s) + 

/i = 0 U=1 L=1 

and let iU(s) = e^^H{s)/T{s), regarded as a function of s e C. It is then shown that 

i/(x) = ^reSsUr (31) 

S 

where res^ denotes the complex residue of F^. at a pole s e C, and the sum is taken over 
all poles s. To ensure convergence, the poles are arranged in an appropriate order [Wri49] 
which is somewhat delicate in general; in our particular setting we will find below that the 
summation is absolutely convergent. 
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To avoid the singularity of g' at zero, let y{x) ^ g\{x + e), so y satisfies 
X~^y{x) — y{x + 1) + y'{x + 1) = 0, —e<x<co. 
This clearly corresponds to (30) with m = n = 1, 


i /=0 J /=1 

fj,=o ( 1/A 0 


M=i \ -1 


b„ = 


M=' 


=0 / 0 


\ 1 


t{s) = 1/A — + 56^. 


(32) 


Note that t'(s) = e^s and r"(s) = e^(s + 1), so if r(s) = 0 then s is a simple root unless 
s = 0, in which case A = 1 and s = 0 is a double root. From (29) we have initial data 


oX + e 


/ \ \ I for a: + e e 10 , 11 , 

y{^) 9 \{x + e) for a: + e e [ 1 , 2 ], 


which we use to evaluate 
’1 


Hr{s) = e 


f'- 


y(u) + y'(u)]e du = 


- e" - e" 


H(s) = Ffi(s) + e 


s+e 


A(l-.) ’ 

where iTi(s)|s=i is understood to be —ee^. If r(s) = 0 then iT(s) = 7 ^ 0, so the poles 

of Fa,(s) = e®^iT(s)/r(s) correspond precisely to the zeroes of t(s). We compute 


e^^H(s) 


:,{x+e)s 


\im{z - s)F^{z) = 
reSs Fx = < , r (s) s 

lim —[z‘^Fx{z)] = 2(x + e + 1/3) 
z^odz^ ^ ^ ^ ’ 


if s is simple pole, 
if A = 1, s = 0 


(recall from above that these are the only two cases for roots of r). Applying (31) gives 


S'JC xs 

9 \{x)= V — for A e ( 0 , 1 ), ^i(a;) = 2 (a; + 1/3) + V —, 


ssC: 

t(s)=0 


ssC\{0}: 

t(s)=0 


(33) 


modulo issues of convergence to be addressed in the next subsection. 


6.3. Roots of r and Lambert W function. The roots of r can be expressed in terms of 
the Lambert W function, which has been very well studied and which we now briefly review.^ 
The number of rooted spanning trees of the complete n-vertex graph is rF‘~^ (Cayley’s 
formula). The tree function T is the associated exponential generating function. 


m = T 


00 IT) _ 1 


71=1 


n\ 


z^. 


From Stirling’s formula (n! '\/27rn(n/e)"'), the sum converges for all \z\ ^ 1/e. It satishes 

the relation T{z) = 2 :exp{T( 2 :)} (see e.g. [Sta99, Proposition 5.3.1]). The Lambert W function 
is dehned by the equation 

z = W{z)e^^^l (34) 

This is a multivalued function, with branches Wk naturally indexed by the integers fc e Z; 
see [CGH+96, §4] and Fig. 3. 


^This function also arose in the analysis [MW13] of the subcritical regime of the minimum mean-weight 
cycle, for seemingly different reasons. It figures prominently in the analysis of random graphs near the phase 
transition, e.g. [JKLP93]; see [CGH+96, §2] for additional applications. Some of the discussion in this section 
is adapted from [MW13]. 
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branch k=2 



branch k=1 


branch k=-1 


principal branch k=0 



branch k=-2 


Figure 3. (This is a dupiicate of [CGH+96, Fig. 4].) The Lambert W function is defined 
as the (muiti-vaiued) inverse of the function f(w) = we^ . The curves in the figure show the 
preimage set /“^(K<o) — they are a subset of M u {rc : ^ = —77 cot 77 } where ^ = Rezc and 
77 = Imre. The curves naturaiiy partition the rc-piane into branches indexed by /c e Z (the 
fc-th branch is the image of Wk)- Foiiowing convention [CGH+96], branches are defined to 
be ciosed in the direction of increasing ImTC. 

The principal branch Wq satishes Wo{z) = —T{—z) for \z\ ^ 1/e, and can be dehned 
elsewhere by analytic continuation. It is straightforward to check that T(l/e) = 1. For d e C 
with \e~^\ ^ 1, we can use T(l/e) = 1 and the relation T{z) = 2 :exp{T( 2 :)} to deduce that 

T{l/e^+^) = l-V^ + 0{6). (35) 

The branch ITo has a cut on 2 ;e(— 00 ,—1/e); and in view of (35), it has a square-root-type 
singularity near ^ = —1/e. Any other branch Wk {k 7 ^ 0) has a cut on z e (— 00 , 0) with a 
logarithmic singularity near z = 0. 

Recalling (32), the solutions of r(s) = 0 are given by 

Sk = l + Wk, Wk = Wk{-l/{eX)), keZ. (36) 

For A e (0,1], Sk is obtained by evaluating Wk precisely on the branch cut. Following the 
convention that branches are closed in the direction of increasing Imre ([CGH+96], see also 
Fig. 3), we see that Sk and s_k-i are complex conjugates for each k ^ 0. When A = 1, Sq and 
s_i are obtained by evaluating Wq and bF_i at the branch point —1/e, giving the double root 
So = s_i = 0. There is a convergent series expansion [CGH^96, eqn. 4.20] for each branch 
Wk] truncating the series and evaluating at —l/(eA) gives, for positive integers /c, 

Sk = - log[7rA(4/c -b l)/2] -f- i7r(4/c -f l)/2 -b 0(1 -b k^Hog k), 

S-k-i = — log[7rA(4/c -b l)/2] — iTi{Ak -b l)/2 -b 0(1 + A;“^log k). 

Thus, for any a: e M there exists a finite constant such that | exp{sfcx}/sfc| ^ Cx/\k\^^^ for 
all k E 7 j, excluding the case sq = s_i = 0. It follows that for a; > 0, the summations in (33) 
are absolutely convergent as claimed. 

6.4. Asymptotics for A near 1. We now extract the asymptotics of gx{x) when A e (0,1] 
is close to 1, and x is large. Recalling that Sk and s_fc-i are complex conjugates for k e 
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Figure 4. Complex roots of t(s) = 1/A — e®(l — s), shown for A = 0.95. 



Figure 5. Taking A = 0.95, the figure shows gosdix) with its one-term and two-term 
approximations in the series expansion (38) (i.e., using the rightmost or two rightmost 
conjugate pairs of Sk from Fig. 4). The three curves are nearly indistinguishable except for 
negative or very small x. 


we can re-express (33) as 

9\ix) = fxix) + 2 2 Re(e^''7sfc), fxix) 

k^l 


j2Re(e"°7so) if A e (0,1), 
[2(0;+ 1/3) ifA = l. 


(38) 


It is strongly suggested by (37) and Fig. 4 that fx{x) is a good approximation to gx{x) in 
the limit of large positive x. We shall prove for A sufficiently near 1 that this is indeed true, 
which amounts to proving that Re Sk is strictly maximized over k e at fc = 0. Recalling 
(36), it is clearly equivalent to prove that Rewk is strictly maximized over k e at /c = 0. 
For A near 1, wq lies near IF (—1/e) = —1. It is clear from the dehnition of the branch cuts 
([CGH'*’96, §4] and Fig. 3) that for A e (0,1] and fc ^ 0, 

\m.Wk e [2/c7r, {2k + l)7r) , (39) 

so \wk\ ^271 for k ^ 1. Combining with (34) gives 


Rewk = — log |eA| — log \wk\ ^ —1 — log(27r) — log |A| ^ —2 for |A| near 1. 

Combining with (37), we see that there is a hnite constant C such that for all A in a 
neighborhood of 1, \gx{x) — fx{x)\ ^ Ce~‘^^ for all a; > 0. 
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We next appr oximate fx{x) for A = e 5 a small positive real. Taking z = —l/(eA) and 
p = V2(l + ez), for \p\ < there is a convergent series expansion [CGH+96, eqn. (4.22)] 

So = 1 + Wo{z) = V pip^ = p-\p^ + ^P^ + 0{p ^). 

Since p = i-\j2{e^ — 1) = ix/^\l + 0(5)], the hrst term in the expansion matches what we 
have already noted in (35). Therefore 

Re So = (25/3) [1 + 0(5)], Imso = V^[l + 0(5)], 

and as a resnlt, for A = with 5 > 0 small, 

/A(a;) = 2exp{(ReSo)ic}sin (Imso)x + arctan[(Reso)/(Imso)] j\so\ 

= 2 exp{(25/3)a;[l + 0(5)]} sin [1 + 0(5)]V^(x + 1/3) j. 

In particnlar, taking 5 | 0 with x hxed, we recover /i(x) = 2(a: + 1/3). 

Lastly we identify the valne of H for which gx, X = e~^, gives the principal eigenfnnction. 
Recalling the discnssion of §6.1, H + 1 corresponds to the smallest positive zero of g\. From 
the above, f\{xo) = 0 for 


xo = 


71 — arctan[(Re So)/ (Im sq)] tt 


3 + 0{VS) . 


Im So 

Recalling that fx approximates gx, we can write 

gx{xo + x) = 0(exp{-f2(xo + x)}) + fx{xo + x) 

= 0(exp{—f2(a:o + 3:)}) + 2exp{(Reso)(a;o + x)} sin[(Imso)a:]. 

From this we see that we can choose x* = Xq + exp{—f2(xo)} and a sufficiently large constant 
x' such that gx has a root at x*, and is non-vanishing between x' and x*. To rule out zeroes 
on [0,x'], note that gi is non-zero on all of M^o? since no hnite H has an eigenvalue of 1. 
It follows by continuity in gx that for A sufficiently near 1, gx is non-vanishing on [0,x']. 
Therefore 


H = Xi, - 1 = 


TT 




3+0(VS). 


5 = 


71 ^ 1 + 

y (i/ +4/3)2 
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